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ABSTRACT 


The  common  assumption  that  the  stress  is  infinite  at  the  tip  of  an  in-plane  crack  is  incon 
sistent  with  the  basic  historical  solutions  for  stress  for  cases  from  which  crack  formulae  have 
been  evolved.  The  latter  formulae  do  not  satisfy  boundary  conditions.  An  appropriate  definitic| 
of  a crack,  as  does  one  presented  here,  should  make  it  obvious  that  such  conditions  are  to  be 
fulfilled  and  when  they  are,  the  meaning  of  stress  intensity  factor  as  the  coefficient  of  a 1/^/r" 
singularity  is  altered.  It  no  longer  represents  an  infinity  of  stress  and  its  connection  with  actual 
failure  stress  through  a stress  concentration  factor  leads  to  a fixed,  rather  than  experimental, 
connection  between  Modes  I and  II  stress  intensity  factors.  Furthur  discussion  of  appropriate 
representations  of  cracks  for  shear  and  normal  stress  loading  is  warranted,  as  well  as  of  tough- 
ness definition  consistent  with  failure  mechanisms  and  with  elastic-plastic  solutions  for  stress. 
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INTRODUCTION 

This  paper*  is  concerned  with  an  apparent  oversight  in  the  usual  current  operational  definitions 
of  stress  intensity  factor  for  in-plane  cracks,  and  in  its  resolution.  This  oversight  leads  to  an  incon- 
sistency between  the  definitions  of  stress  intensity  factor  for  Modes  1 and  lit  as  well  as  to  a lack  of 
definition  for  stress  intensity  factor  for  combinations  of  these  modes,  assuming  the  intent  of  the 
definitions  is  to  make  fracture  stress  a unique  common  basis.  The  oversight  stems  from  a lack  of 
satisfaction  of  crack  boundary  conditions.  This,  in  turn,  derives  from  setting  the  minor  axis  of  the 
Inglis  solution  equal  to  zero,  in  evolving  the  formulation  of  the  current  definitions  of  stress  intensity 
factor.  The  consequent  aberration  docs  not  occur  if,  in  reducing  the  Inglis  solution,  the  tip  radius  is 
retained  along  with  the  crack  length,  even  though  the  radius  may  be  assigned  as  small  (but  finite)  a 
value  as  one  chooses.  There  is  no  infinity  of  stress. 

Thus  the  radius  becomes  a part  of  our  suggested  definition  of  a crack  and  accordingly  cracks 
can  be  regarded  as  notches. 

Before  discussing  the  proper  deduction  of  stresses  near  crack  tips,  we  will  briefly  review  the 
current  operational  definitions  of  stress  intensity  factor  and  show  that  they  are  not  inconsistent 
with  a crack  tip  radius  concept  and,  in  fact,  can  even  be  regarded  as  a stress,  not  the  usual  stress 
X y'Tengtfi,  so  long  as  fracture  stress  failure  of  the  cracked  material  is  envisioned. 

This  paper  is  an  elaboration  of  part  of  a former  one  [ 1 1 in  which  the  writer  gave  his  own  defi- 
nition of  a crack  and  indicated  how  he  had  been  treating  cracks  as  notches  in  the  prediction  of  failure. 

On  Current  Definitions  of  K 

According  to  current  complex  variable  analyses  which  assume  rectilinear  elasticity,  the  stress 
intensity  factors,  K,  for  Modes  I and  II  are  determined  as  the  coefficient  of  a l/y/2r  singularity  by 
the  formula 

Kj-iKn  = 7y/l  lim  Z1/2  (Z)  Z-  0 (1) 

in  which  <p  is  a stress  function  whose  expansion  commences  with  a-^/Z  term,  r is  the  radial  distance 
from  the  crack  tip. 

We  now  consider  Mode  I,  i.e.,  crack  axis  is  perpendicular  to  the  direction  of  the  resultant 
loading. 

The  above  definition  really  stems  from  taking  K to  be  the  coefficient  of  (2r)'1/2  f(0)  in  the 
expressions  for  the  rectangular  components  of  stress  as  distance,  r,  from  the  tip  of  the  crack  ap- 
proaches zero,  i.e.,  typically^ 


(2) 


•Manuscript  received  3 June  1975. 

fMode  I:  Loading  stresses  at  infinity  perpendicular  to  and  along  the  crack  axis;  Mode  II:  Loading  stresses  at  infinity  are  shearing 
stresses  along  and  perpendicular  to  the  crack  axis. 

^ The  complete  expressions  may  be  found  on  pages  7 and  1 0. 


I.  BEEUWKKS,  R.,  Jr.  Analysis  o/  Failures.  Proceedings  of  the  Third  Sagamore  Ordnance  Materials  Research  Conference, 

December  1956.  co-sponsorcd  by  the  Otdnancc  Materials  Research  Office  (now  Army  Materials  and  Mechanics  Research  Center) 
and  the  Office  of  Ordnance  Research  of  the  U.  S.  Army,  p.  99. 
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which  to  be  dimensionally  correct,  and  because  stress  is  proportional  to  load,  may  also  be  written* 


i.e., 


(k/2)Sn  \/  a 

S = f (A) 

s/2r 


K = (k/2)  Snv/i 


where  Sn  is  a nominal  loading  stressf , “a”  is  a length  which  directly  or  indirectly  characterizes  crack 
size  and  k is  a function  of  crack  specimen  geometry.  Thus  for  the  standard  reference  case,  a crack  in  a 
plate  many  crack  lengths  long  and  wide,  and  loaded  cross-wise  to  the  crack  by  a uniform  end  loading 
stress  S„,  k = 2 if  “a”  is  half  the  crack  length.  If  this  plate  were  reduced  symmetrically  to  width  w, 
keeping  the  crack  in  the  center,  k as  2 (cos  ira/w)-l/2.  If  the  original  plate  had  been  sectioned  along 
the  specimen  axis,  cutting  the  crack  in  half,  and  the  meaning  of  Sn  .and  “a”  were  retained, 
k = (2X1 .1 2147).  If  instead  of  the  center  crack  our  original  specimen  had  two  very  (~  infinitely) 
deep  cracks  symmetrically  situated  to  each  other  on  opposite  sides  of,  and  perpendicular  to  the 
specimen  axis  with  the  nominal  stress  conveniently  selected  to  be  the  axial  load  divided  by  the  net 
section  area  between  the  crack  tips  and  “a”  selected  to  be  half  the  distance  between  the  tips,  then 
k = 4/ir. 

This  formula  for  stress  purports  to  uniquely  depict  the  stress  field  immediately  about  the  crack 
tip;  k embodies  the  effect  of  specimen  geometry  and  load  distribution  but  not  the  distribution  of 
stress  at  the  tip. 

Clearly  the  formula  may  also  be  written,  by  dividing  both  numerator  and  denominator  by  s/p. 


S = KV?_.  f(fl)  = — 

s/Trfp  y/2  r/p 


in  which  we  see  that  distance  is  now  measured  in  terms  of  a 
radius,  i.e.,  distance  is  now  measured  by  r/p. 

Also 


f(0)  (3) 

characteristic  length,  such  as  a crack  tip 


K0  = K Vp  = (k/2)  SnV/Sfc 


(4) 


a quantity  which  has  the  dimension  of  a stress  and  an  appearance  like  that  of  the  concentrated  stress 
at  the  base  of  a notch.  It  could  be  computed  if  p were  a known  quantity. 


•k/2  i«  written  u \/\/n  in  tome  literature  (Y  depends  on  the  dimension  chosen  for  "a").  The  reason  for  choosing  k/2  for  the  constant, 
instead  of  simply  k,  is  to  identify  k with  the  stress  concentration  factor  coefficient  of y£/p  as  will  be  seen  later  on.  Also,  our  K\^t~ 
is  often  called  K,  as  a result  of  an  energy  approach,  but  is  redundant  and  distracting  in  the  formulae  for  stresses, 
tin  determining  K for  a material,  Sn  and  “a”  are  measured  at  a failure  load  defined  by  some  Standard,  especially  the  ASTM  Standard 
for  toughness  testing. 


2 


Now  the  term  “stress  intensity  factor”  when  used  to  evaluate  toughness  of  materials  may  be 
taken  to  imply  that  a fracture  stress,  S = F,  is  attained  at  failure  according  to  the  formula  for  K in- 
volving stress  and  crack  length  evaluated  at  failure.  That  is,  since  f (0)  = 1 when  S is  a maximum. 


S = F = K0A/2r7p 


(5) 


where,  to  have  similarity  at  failure  in  different  tests,  r/p  = constant.  If,  for  example,  r/p  - 1/2, 

F = KVp  = K0 

or,  if  F and  K are  known  experimentally 

p = (K/F)2 


Instead  of  K,  G is  sometimes  used,  where,  in  our  first  example  referring  to  a central  crack  in  a 
tension  member 


S 


n 


GE 

7ra(l-p2) 


(6) 


and  G corresponds  to  the  failure  value  of  Sn  when  the  crack  length  is  2a.  Here,  in  analogy  to  the 
Griffith  formula. 


S 


n 


2 E-y 
7r  a(  1-p2) 


for  crack  failure  when  surface  energy  is  controlling,  writers  assume  G represents  the  energy  absorbed 
at  crack  failure,  per  unit  length  of  crack  growth. 


This  assumption  is  unwarranted  if  G and  K are  computed  from  the  same  critical  values  of  Sn 
and  “a”.  While  an  energy  criterion  must  be  met,  it  is  not  necessarily  controlling.  However,  assum- 
ing G is  controlling,  we  have 


and 


= k2  (7) 

jrl  1-p2) 


_ = /K\2  _ EG 

P \fJ  tt(1-p2)F2 
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Definition  of  a Crack  | 1 1 

In  accordance  with  our  previous  remarks,  we  assume  that  for  the  mathematical  treatment  of 
macroscopic  cracks,  the  crack  may  be  regarded  as  a symmetrical  notch  with  a smoothly  varying 
curvature  and  a parabolic  tip  of  sufficiently  small  curvature  compared  with  its  depth.  “Sufficiently 
small  curvature”  means  that  an  elasticity  stress  analysis  would  show  that  the  peak  stress  in  the  re- 
gion of  the  tip,  is  inversely  proportional  to  the  square  root  of  the  tip  radius  of  curvature  to  any 
close  approximation  desired.  Since  any  curve  with  continuous  second  derivatives  is  characterized 
at  a point  by  its  curvature,  it  will  be  seen  that  this  definition  includes  more  than  obviously  crack- 
like ellipses  and  hyperbolae. 

This  definition  clearly  requires  satisfaction  of  boundary  conditions  at  the  tip  of  the  crack. 

This  is  not  true  of  the  usual  conception  of  a line  crack  (at  least  in  the  way  it  has  been  treated).  In 
the  latter  case  any  symmetrical  se'f-equilibrating  combination  of  concentrated  loads  could  be  ap- 
plied to  the  tip  along  the  crack  ax  >,  and  satisfy  the  boundary  conditions  as  well  as  does  the  analysis 
one  has  come  to  expect,  which  m;  kes  the  stress  in  the  crack  axis  direction  equal  to  that  perpendicu- 
lar to  the  crack  axis,  instead  of  ze  o*  as  it  should  be  and  would  be  in  the  case  of  the  parabolically 
tipped  crack.  These  loads  would  also  cause  stress  to  decrease  with  distance  from  the  crack  tip  at  a 
different  rate  than  the  inverse  square  root  ol  this  distance,  in  contrast  to  the  current  treatments. 

This  definition  excludes,  or  relegates  to  separate  treatment,  a crack  with  parallel  sides  and 
circular  tip.  It  and  other  shapes  should  be  studied  separately  not  only  for  pedagogical  reasons,  but 
because  they  may  conceivably  be  encountered,  or  prove  useful  experimentally.  With  them  peak 
stresses  are  reached  at  different  loads  and  amount  of  yielding  near  the  crack  tips.  However,  the 
definition  adopted  here  derives  historically,  as  do  the  conventional  crack  formulae  for  K,  from  the 
well-known  Inglis  classical  solution  for  stresses  about  elliptical  holes  and  thus  may  be  readily  visual- 
ized and  the  Inglis  solution  utilized  by  the  reader  in  deriving  the  results  to  be  presented  below.  The 
definition  is  however,  not  only  suitable  for  present  purposes,  but  separate  investigation  [2] , not 
presented  here,  indicates  that  its  use  is  not  limited  to  smooth  contours  but  that  the  material  to 
which  it  is  applied  may  be  considered  to  have  an  effective  radius  when  there  is  some  plasticity  about 
the  crack  tip,  an  effective  radius  which  is  a material  property  and  is  associated  with  failure  below 
the  crack  tip  with  its  irregular  contour. 

Stress  Formulae  with  Parabolic  Crack  Tip 

The  general  formulae  for  stress  components,  in  the  vicinity  of  the  crack  tip  may,  as  previously 
indicated,  conveniently  be  deduced  from  the  Inglis  solution  for  stresses  about  an  elliptical  crack,  if 
this  solution  is  expressed  in  terms  of  the  semi-major  axis  “a”  and  the  radius  of  curvature  at  its  tip.  p, 
instead  of  “a”  and  the  semi-minor  axis  “b”.  That  is,  for  “b”  we  first  substitute  b , then  re- 

place the  coordinate  system  by  a natural  one  that  is  near  the  tip  and  finally  reduce  the  resulting 
formulae  to  the  limiting  case  where  a/p  >>  1 and  distance  from  the  tip,  measured  in  tip  radii,  is 
small. 

The  resulting  formulae  for  the  rectangular  components  arc  as  follows.  The  corresponding 
formulae  for  the  polar  and  parabolic  coordinate  components  are  submitted  in  the  Appendix.  The 
Airy  stress  function  is 

* Assuming  no  internal  pressure  in  the  crack 

2.  BHUWKIS,  R.,  Jr.  Characteristics  of  Crack  Failure.  Surfaces  and  Interfaces,  Syracuse  University  Press  Syracuse  N Y 
v.  II,  1968,  p.  277. 
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p,  - <k, s,/2)  i.mismi  +£2di^i  l8| 

(2  r/a )3/2  (2  r/ap/2 

Because  of  the  p term,  the  stresses  do  not  decrease  witli  distance,  r,  from  the  tip  as  r -1/2. 

The  origin  of  rectangular  coordinates  (x,y)  and  polar  coordinates  (r,  0)  is  at  a distance  p/2  be- 
hind the  nose  of  the  crack-like  notch,  on  the  x axis,  p is  the  radius  of  curvature  of  the  tip  of  the 
nose,  x is  along  the  crack  axis  the  negative  direction  lying  completely  within  the  crack.  The  posi- 
tive y axis  is  counter-clockwise  from  the  positive  x axis.  The  polar  coordinate  angle,  0,  is  positive 
in  the  direction  from  positive  x to  positive  y.  The  coordinates  (x/a),  (y/a)  and  (r/a)  are  oriented  and 
centered  in  the  same  way  as  (x,  y)  and  r. 

The  equation  of  the  parabola  is,  in  polar  coordinates, 

(r/p)  (2  cos2  0/2)  = (r/p)  ( 1+  cos  0)  = 1 (9) 

and  in  rectangular  coordinates, 

x/p  = (1/2)  ( 1 - (y/p)2] 

from  which  it  is,  indeed,  evident  that  the  tip  is  at  a distance  of  x = p/2  in  front  of  the  origin.  It  is 
also  obvious  that  with  distance  measured  in  tip  radii,  all  parabolae  have  the  same  size  and  appear 
crack-like  if  x > > p. 

Mode  I:  “Normal  Stress”  - Crack  Tip  Stress  Symmetric  with  Respect  to  Crack  Axis 

Only  the  loading  stress  perpendicular  to  the  axis  appears  in  the  expressions  for  stress  near  the 
crack  tip  if  a/p  » 1 , as  assumed  here. 


2Sx 

(p/a)  cos  3 0/2 

cos  0/2(1 -sin  0/2  : 

k|S| 

(2  r/ap/2 

(2r/a)  1/2 

2 sy  _ 

(p/a)  cos  3 0/2 

cos  0/ 2(  1+  sin  0/2 

k]  S! 

(2  r/ap/2 

(2  r/a)  1/2 

2Sx*  _ 

(p/a)  sin  30/2  + 

cos  0/2  (sin  0/2  cos 

kj  S, 

(2  r/a  p/2 

(2r/a)l/2 

The  stresses  on  the  boundary,  p/a  = (2  r/a)  cos-  0/2,  are,  keeping  the  contributions  of  each 
term  separate  before  adding  them  up, 
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- (cos^  0/2  cos  30/2)  y/a/p  + cos2  0/2  ( I -sin  0/2  sin  30/2)  \A/P 
+ (2  si n ^ 0/2  cos2  0/2)  y'a/p 

+ (cos3  0/2  cos  30/2)  ^ajp  + cos^  0/2  ( l'+sin  0/2  sin  30/2)v/a/p 
(2  cos4  0/2)  v/a7*r 

- (cos3  0/2  sin  3 0/2)y/j]p  + cos2  0/2  (sin  0/2  cos  30/2)^afp 

(-2  sin  0/2  cos3  0/2)  y/a/p 

However,  the  stresses*  in  the  directions  u,  v of  parabolic  coordinates. 

u2  = (2  r/a)  cos2  0/2,  v2  = (2  r/a)  sin2  0/2,  boundary  p/a  = (2  r/a)  cos2  0/2 

are  more  useful  in  determining  the  relative  contributions  to  the  boundary  stress  state. 
These  are 

2 Su  _ -p/a  cos  0/2  + cos3  0/2 
kI  Si  (2  r/a)3/2  (2  r/a)1/2 

so  that  on  the  boundary 


2bx 

kfSI 


= 


ki  Sj 


2 ^xy 

k,S| 


r-lr-  = (-cos'*  0/2)  y/ajp  + (cos'*  0/2)\/aJp  = 0 

Kpj 

2 Sv  p/a  cos  0/2  cos  0/2  (2  - cos2  0/2) 

kisT  = 12^5372"  72  r7a)l75— 


(ID 


(12) 


03) 


* The  following  formulae  are  conveniently  derived  using  parabolic  coordinates.  However,  they  may 
be  found  from  the  rectangular  components  of  stress  by  taking  account  of  the  fact  that  in  an  ele- 
mentary right  triangle  with  hypotenuse  ds  along  the  parabolic  boundary  the  angle  between  dy 
and  ds  is  0/2. 
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so  that  on  the  boundary 


2 S 

^7=  (cos4  0/2)  Va/p  + (2  cos2  0/2  - cos4  0/2)  y/JJp 
= (2  cos2  0/2)  v/a7p 

2 Suv  _ -p/a  sin  0/2  cos2  0/2  sin  0/2 

kISI  (2r/a)3/2  (2r/a)l/2 

so  that  on  the  boundary 

2 Syy  ^ _ 

= (-cos-3  0/2  sin  0/2)  ^afp  + (cosJ  0/2  sin  0/2)  y/ajp  = 0 


Here  again,  the  contributions  of  the  two  types  of  terms  have  been  kept  separate  before  adding 
them  up. 

In  expressions  for  Sy  and  Sv  for  0_=  0,  we  recognize  the  usual  formula  for  stress  concentration 
of  a very  sharp  elliptical  notch,  2^/a/p,  where  kj  = 2 as  it  is  for  that  case  when  Si  is  the  loading 
stress  far  from  the  notch. 

As  stated  previously,  the  first  term  was  unfortunately  omitted  in  the  definition  of  stress  inten- 
sity factor  so  that  the  formulae  for  stresses  were  written 


cos  0/2  (1-sin  0/2  sin  3 0/2)  (14) 

727)1/2 

K]  cos  0/2  (1+  sin  0/2  sin  3 0/2) 

(27)7/2 

K|  cos  0/2  (sin  0/2  cos  3 0/2) 

(27)172 

Thus  if  these  latter  formulae  were  evaluated  on  the  boundary,  at  0 = 0 where  the  stress  is  maxi- 
mum, the  term  retained  in  them  contributes  but  half  the  true  stress  which  is  given  by  the  prior  for- 
mulae. However,  comparing  these  sets  of  formulae,  it  is  evident  that 


ki  Si 

~2“  \Z'd  ~ Kl  (15) 


1 


and  hence  the  expression  for  K|  computed  mathematically  as  the  coefficient  of  (2  r)-*/2  can  bc 
used  to  compute  the  coefficient  of  stress  concentration.  k|,  by  tins  relation,  without  error,  even 
though  an  improper  expression  for  stress  has  been  used  in  defining  K. 

If  S is  assumed  to  be  F,  where  F is  the  “nil-ductility”  tensile  fracture  stress  and  K = K|c 


F=k[S|  s/u/p  = 2 K| e/x/P  06) 

by  the  above  relation,  i.e., 

P = (2  K,c/F)2 

This  p could  be  a machined  in  radius  and  the  stress  distribution  might  be  rectilinearly  elastic 
up  to  S = F,  so  that  this  case  would  be  one  of  classic  brittle  failure.  On  the  other  hand,  this  p might 
be  a unique  characteristic  length  obtained  by  this  formula  from  any  crack  Kjc  test,  a length  which 
is  in  some  cases  unequivocally  a crack  tip  radius. 

Mode  II:  "Shear”  Crack  Tip  Stress  Anti-Symmetric  With  Respect  to  Crack  Axis 

The  resultant  component  of  the  shear  loading  which  is  parallel  to  the  crack  axis,  is  here  re- 
garded as  a clockwise  couple  whose  forces  lie  symmetrically  above  and  below  the  crack  axis,  so  that 
the  force  above  points  away  from  the  tip  in  the  + x direction. 

The  Airy  stress  function  is 

Fn  = (kn  Sri)  Li£/i_5!!LW?>  + 3 sin  0/2  cos-  0/2") 

" " n)  [ <2r/„)3/2  (2  r/a)l/2  ll7’ 


The  stresses  are 


Sx 

+ (p/a)  sin  30/2)  _ 

sin  0/2  (2  + cos  0/2  cos  3 0/2) 

(18) 

kII  Sji 

(2  r/a)3/2 

(2  r/a) I/2 

sy  = 

_ (p/a)  sin  3 0/2  + 

sin  0/2  (cos  0/2  cos  3 0/2) 

kll  Sn 

(2  r/a)3/2 

(2  rOj T^~ 

^xy 

(p/a)  cos  3 0/2 

. cos  0/2  (1  - sin  0/2  sin  3 0/2) 

kll  SII 

(2  r/a)  3/2 

72  r/a)*/2 

As  was  the  case  with  Mode  1,  the  stresses  in  the  directions  u,  v of  parabolic  coordinates  are 
more  useful  in  determining  the  relative  contributions  to  the  boundary  stress  state. 
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I 


These  are 


= 0>/a)  sin  0/2  sin  0/2  cos2  8/ 2 

kII  SII  (2r/a)3/2  (2  r/a)>/2 


(19) 


so  that  on  the  boundary 


Su  _ cos3  0/2  sin  0/ 2 \JJJp  - sin  0/2  cos3  6/2  y/'djp  =0 

Mi  sn 


^ = _ (p/a)  sin  812  _ sin  8/2  (1  + sin?  8/2) 

k|j  S||  (2  r/a)3/?  (2r/a))/2 


so  that  on  the  boundary 


i , 

fjTsjT cos:*  s'n  v^a/p  - (2  cos  0/2  sin  0/2  - cos3  0/2  sin  8/2)\[dJp 

_ _ s»n  812  cos  2 0/2  _ sin  0 (2  - cos2  0/2)  y/y^- 


sin  0y/i7p 


Suv 

Ml  S|I 


(p/a)  cos  0/2  + cos3  0/2 


(2  r/a)3/?  (2  r/a)>/2 

so  that  on  the  boundary 
S 

j-~-  = “ cos4  0/2  ^a/p  + cos4  0/2  y/a/p”' 


= 0 


i 


p. 


As  before,  the  contributions  of  the  two  types  of  terms  have  been  kept  separate  before  adding 
them  up. 

In  the  expression  for  Sy  we  recognize  the  formula  for  stress  concentration  of  a very  sharp 
elliptical  notch  *y/a[p  , corresponding  to  extremes  of  sin  0,  0 = ± 7r/2,  where  k|j  = I as  it  is  for 
the  case  when  S|j  is  the  loading  stress  far  from  the  notch. 

Again,  as  we  observed  previously,  the  first  term  was  unfortunately  omitted  in  the  definition  of 
stress  intensity  factor  so  that  the  formulae  for  stresses  were  written 


K||  sin  0/2  (2  + cos  0/2  cos  3 0/2) 


(2  r)*/2 


Sy  - + 


K 1 1 sin  0/2  (cos  0/2  cos  3 0/2) 

(2  r)>/2 


^xy  + 


K|j  cos  0/2  (1  - sin  0/2  sin  3 0/2) 
‘(2r)T/2 


Thus  if  the  latter  formula  were  evaluated  for  stress  along  the  boundary  p = 2 r cos2  0/2,  we 
would  have 


Sy  = Sy  cos2  0/2  + Sx  sin2  0/2  - 2 Sxv  sin  0/2  cos  0/2 


Sfc|V^  =~  { ~ cos3  °/2  cos  3#/2)  + 2 sin2  0/2 


+ sin2  0/2  cos  0/2  cos  3 0/2) 

+ 2 cos2  0/2  - 2 sin  0/2  cos2  0/2  sin  3 0/2)| 

— - | 2 - 2 cos3  0/2  cos  3 0/2  + cos  0/2  cos  3 0/2) 

- 2 sin  0/2  sin  3 0/2  cos2  0/2  | 
syiJj_  1 2 + cos  0/2  cos  3 0/2  - 2 cos2  0/2  cos  0 j 

S*2 ""  | 2 + cos  0/2  cos3  0/2  - 3 cos  0/2-4  cos3  0/2+2  cos  0/ij| 

sin  0 ( - 1 1 sin  0 ( , „) 

- , < 2 - cos-  0/2  > = “4  \ 3 - cos  ® | 
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Hero  we  have  the  same  expression  that  appeared  above  as  the  coefficient  of>/a Jp  in  the 
second  term  ot  Sv  in  parabolic  coordinates.  This  expression  has  a maximum  at  0 = n/2  + cos- 1 
1(3  - V 1 7i/41  but  this  is  misleading.  The  true  extremes  are  at  + 0 = n/2  as  was  seen  in  the  total 
expression  for  Sv  Thus  the  above  expression  is  + 3/4  at  0 = ± tt/2,  instead  of  + 1 , the  value  of  the 
complete  expression. 

Unfortunately  for  consistency  of  K in  representing  stress,  this  differs  from  the  ratio  of  1 to  2 
found  for  Mode  I. 

For  Mode  II,  it  appears,  in  comparing  the  above  set  of  formulae  expressed  in  terms  of  Kii  with 
that  satisfying  the  boundary  conditions,  that 


kji  Sji  * Kjj  (22) 


and  hence  the  expression  for  Kjj  computed  mathematically  as  the  coefficient  of  (2  r)-1/2  can  be 
used  to  compute  the  coefficient  of  stress  concentration,  kjj,  by  this  relation,  without  error,  even 
though  an  improper  expression  for  stress  has  been  used  in  defining  K,  an  expression  which  does  not 
even  specify  the  angle  6 at  which  maximum  stress  occurs. 

Again,  if  S is  assumed  to  be  F,  where  F is  the  nil  ductility  tensile  fracture  stress  and  K|j  = K|jc’ 


F - kji  Sjj  v/a/p  = kjic/ y/ p (23) 


by  the  above  relationship  of  kjj  Sjj  y/H  to  K|| 

Therefore 

p = (K,Ic/F)2 


As  stated  with  reference  to  Mode  I,  thus  pcould  be  a tnachined-in  radius  and  the  stress  distribu- 
tion might  be  rectilinearly  elastic  up  to  S = F so  that  this  case  would  be  one  of  classic  brittle  failure. 
On  the  other  hand  thisp  might  be  a unique  characteristic  length  computed  by  this  formula  from 
any  crack  Kj|c  test,  a length  which  is  in  some  cases  unequivocally  a crack  tip  radius. 


Combination  of  Modes  l and  II  : Maximum  Boundary  Stress 

We  take  axes  and  loading  stresses,  Sj  and  Sj|  as  discussed  under  Modes  I and  II  Assume 
Sj  > 0,  Sj|  > 0. 

Then,  from  the  results  already  found  for  the  boundary  stress,  Su  = Suv  = 0 and 


Sv  fcosfl) — ^ ~ (sinfl)  k | j Sjiy/a/p 


I I 


(24) 


which  is  clearly  an  extreme  when 


sin  0/  cos  0 = -kjj  Sjj/(kj  Sj/2),  thus  — ir/2  < 0 < 0 


Thus,  substituting  for  sin  0, 


( kt  Si  (k[j  S || )2  j 

Sv.  ext  = ^ (I  + cos  0 ) --  + — — cos  0j  y/afp 


and  noting  from  the  extremum  condition  that 


sin-  0 ( = I - cos-  0)  = [ kjj  Sjj/  (kj  Sj/2)]2  cos-  0 


so  that 


we  have 


cos  0 = 


k|  Sj/2 


\Z~  (ki  Sj/2)2  + (k|i  Sjj)2 


Vext  ~ { (ki  si/2)  + \/(kj  Sj/2)2  + (kn  Sn)2  l y/ajp 
This  clearly  reduces  to  k|  Sj  \fa[p , 0 if  S||  is  zero  and  + kjj  Sjj  y/a. Tp  if  Sj  is  zero,  as  it  should. 


ki  Si  x/if  _ 

Ki  = 2 ; Kn  = k|i  Sn  y/a 


it  may  also  be  written 


Sv,ext  yfp  ~ 


*i  + VK } + K n 


We  would  certainly  expect  a sharply  notched,  a/p  > > 1,  isotropically  brittle  elastic  material 
to  fail  at  the  notch  surface  when  either  of  these  expressions  reach  a fracture  stress  Sv  = F. 

Since  p = 2 r cos20/2  = r ( I + cos0),  the  radial  distance  to  the  point  on  the  parabola  where 
this  extreme  occurs  is 

r _ 1 

p 1 + coslF  ... 


I + (k|S|/2)/>v/^  (kjS[/2)2  + (knS,,)2 
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* 


■ •« 

I 


i 

ft 


measured  in  tip  radii.  From  the  above  expressions  for  K],  and  K||,  this  may  also  be  written 


or,  if  kj  = 2 k||,  the  normal  case. 


r^ 

P 


S,/V?  + S2 


(28) 


(29) 


Unfortunately,  in  general,  we  cannot  define  or  specify  a mixed  mode  stress  intensity  factor 
from  the  preceding.  This  is  because  the  formulas  for  K in  terms  of  Sv  and  p are  inconsistent  for 
Modes  I and  II.  That  is,  for  given  Sv  and  p,  Kj  and  Kjj  should  be  the  same  since  both  refer  to  a 
stress,  Sv,  acting  at,  and  tangentially  to  the  boundary  of  the  notch,  a stress  which  would  cause  fail- 
ure if  it  reached  a unique  fracture  value,  F,  whether  or  not  it  was  caused  by  Mode  I or  Mode  II  load- 
ing even  though  these  cause  the  maximum  stress  (i.e„  Sv  = F)  to  be  reached  in  different  locations. 

Thus 


Sv 

ki  Si  ^/ajp  = 2 K \ly/p 

Sv  = 

kII  Sn  y/afp  = K n/>/p~ 

so  that  for  the  same  Sv  and  p, 

K,  = 

Syy/p/2 

Kll  - 

Sy  \/p 

which  K are  not  equal. 

The  notch  stress  factors  kS,  i.e.,  kj  Sj  and  k|j  Sjj,  do  not  suffer  from  this  defect, 
value  of  Sy  and  a fixed  value  of  s/a/p 

For  a fixed 

kiSi  = 

kII  Sll 

(30) 

Thus 

kS  = { (k!  S,/2)  +V(kiS!/2)2  + (kn  Sn)2  l 

(31) 

is  the  stress  factor  representing  the  combined  loading  Sj,  Sjj. 
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One  may  well  wonder  why  K instead  of  k was  introduced,  or  at  least  why  the  proportionality 
factor  of  K to  S was  not  made  the  same  for  Modes  I and  II. 

Normally,  with  stresses  S|  and  S||  specified  as  discussed  previously. 


kI  = 2k,i 


so  that  in  this  case 


J-  {si  sr  + sj-  } v^ 

kH  {SI  +\Af  + } v/i/p 


and  we  can  define  a mixed-mode  stress  intensity  factor,  Kj  jj  in  terms  of  an  effective  nominal 
stress,  S|  jj,  i.c., 


sl.  11  = 


sj  +>y/s; 


I SII 


Sv,  max  = — S|  jj  v/a/p  = knSj  jjv/a7p 


SI,II  >A  = kII  SI,II  Va 


KI,  II 


Mixed  Mode  Example:  Tension  Member  with  Central  Elliptical  Crack  at  an 
Here,  if  T is  the  nominal  tensile  stress. 


S|  = T cos^  0 


S||  = T sin  0 cos  0 


where  0 is  measured  counterclockwise  from  the  lateral  axis  perpendicular  to  the  tensile  axis. 


N 


t 


so  that 


Thus 


I - 2; 

k„  = 1 

K1,I1 

= V* 

N1} 

= y/T  j 

( Si  + yf 

= Ty/a 

| cos2  0 

= wr  { 

cos2  0 

^1, II, max  ~ 2 T y'a 

+ sin^  0 


= *1 


Similarly, 


= S1  * ^ 


$2  + «2 

i 


V^Tp 


(36) 


■?  a 

COS'"  0 


(37) 


(38) 


(39) 


(40) 


K,  that  they  are  propor- 
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(k|/2)  S|  ^/a  = K[ 


(41) 


k||  Snv/a  = (k|/2)  Sj|  ,/a  = Kn 


so  that  the  stress  concentration  factors  for  these  two  cases  would  be  the  same  if  p were  the  same, 
i.e..  we  would  have 


(kj/2)  vVp  = k||  y/iTp 


(42) 


However,  the  standpoint  of  this  paper  is  that  this  is  an  improper  approach,  that  the  proper  one 
requires  that  K|  = Kjj  if  the  fracture  stress  is  fixed. 

Or  it  mi^ht  be  objected  that  though  the  p term  has  the  same  dimensions  as  the  (2  r )-'/2  term, 
i.e..  p/( 2 rp/-  vs.  (2  r)-*/2,  that  it  nonetheless  does  diminish  at  a much  greater  rate  than  does  the 
(2  r)-'/-  term,  and  therefore  might  be  regarded  as  not  very  significant,  especially  in  the  non- 
homogeneous  crystalline  state  of  a crack  front  of  metallic  material.  From  this  standpoint  one 
might  also  question  the  (2  r)“H  - term  also,  and  thus  this  whole  elasticity  approach  to  the  problem. 
Where  does  one  decide  what  to  retain?  Such  an  argument,  if  rigorously  adhered  to,  would  ultimate- 
ly question  the  usefulness  of  the  macroscopic  theory  of  elasticity,  as  applied  to  metallic  materials, 
even  in  regions  away  from  serious  stress  concentration.  Thus  in  this  connection  we  point  out  that 
even  a crack  front  ordinarily  passes  many  grains  and  hence  a theory  assuming  homogeneity  may  be 
expected  to  have  some  practical  use.  The  irregular  front  may  even  contribute  to  the  concept  rather 
than  negating  it;  for  example,  there  may  be  an  effective  radius  due,  among  other  possibilities,  to  the 
differences  in  crack  length  and  therefore  openings  from  point  to  point  along  the  front.  In  fact,  as 
indicated  previously,  investigation  |2]  [3]  indicates  that  crack  failure  is  generally  sub-surface  and 
that  an  effective  radius,  such  as  a radius  consistent  with  the  crack  opening  contour  near  the  tip,  is 
suitable  and  sufficient  for  a radius  term  in  the  calculation  of  the  fracture  stress  causing  failure. 

However,  the  objection  of  the  preceding  paragraph  is  not  germane  to  the  main  purpose  of  this 
paper,  though  it  is  to  whether  the  elasticity  treatments  of  K are  useful  or  not.  More  on  this  below. 
The  point  here  is  the  validity  and  self-consistency  of  the  elasticity  treatment  itself,  whether  useful 
or  not.  That  is,  elasticity  theory,  only,  is  used  in  the  formula  for  K which  is  used  in  Kc  determina- 
tions. That  theory  should  be  cast  into  a form  which  satisfies  boundary  conditions.  And  it  should 
embrace  artificial,  machined  in  cracks  of  indefinitely  small,  but  assignable,  tip  radius.  These  factors 
are  accounted  for  in  the  treatment  in  this  paper,  but  it  is  well  to  remember  also  that  Neuber  |4|  in 
his  treatment  of  sharp  notches,  using  a texture  particle  zone  below  the  notch  root  characteristic  of 
the  material  employed,  took  care  to  satisfy  boundary  conditions.  His  result  was  that  a material 
constant,  one  half  the  zone  depth,  i.e.,  e/2  = p ”,  took  the  place  of  the  tip  radius  p in  the  usual 
theory. 

For  self  consistency  the  viewpoints  expressed  above  and  in  the  text  for  crack-like  notches 
must  hold  for  V-shaped  notches,  or  notches  whose  tips  are  V-shaped.  The  crack-like  notches  should 


3.  TAGGART,  R.,  and  WAIII,  K.,  University  of  Washington,  Seattle.  Crack  Opening  Displacement  During  the  Tensile  Loading  of 
Otctile  and  brittle  Notched  Plates.  Private  Communication,  April  1971. 

4.  NfcUBI  R,  II.  Theory  of  Notch  Stresses  Principles  for  Pxact  Stress  Calculation.  The  David  W.  Taylor  Model  Basin,  Translation  74, 
November  1945,  p.  160,  and  reference  to  “Pine  neue  clastichc  Materialkonstantc"  by  L.  Foppl,  Ing -Archiv  v 7 no  4 1936 

p.  229-736.  ’ ’ ’ 
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be  considered  as  a particular  case  of  the  V-shaped  notches.  In  fact,  the  sides  of  the  V should  be  con- 
sidered to  be  the  asymptotes  of  the  actual  notch  shape,  not  the  notch  itself.  Just  as  in  the  crack 
case,  in  order  to  satisfy  boundary  conditions  one  would  incorporate  a tip  radius  into  the  notch  and 
would  then  find  that  no  matter  how  small  the  radius  is  taken  to  be,  that  the  boundary  stress  term 
in  which  it  is  incorporated,  contains  a fixed  proportion  of  the  total  boundary  stress.  One  should 
not  therefore,  take  the  pure  V limit  as  the  V notch  anymore  than  the  particular  case  of  V-notch.  the 
straight  line,  as  the  crack.  From  this  point  of  view,  certain  well  known  solutions  15)  |6]  for  open 
notches  are  incomplete,  and  require  an  additional  term  incorporating  a radius,  to  complete  them. 

That  these  comments  about  tip  radii  are  not  unimportant  will  be  appreciated  when  it  is  con- 
sidered that  one  could  not  solve  for  the  stresses  about  a hole  having  V-shaped  corners  through  the 
use  of  a stress  function  (with  undetermined  coefficients)  which,  while  assuming  the  corner  is  reached 
as  the  zero  limit  of  a radius,  assumes  that  the  corner  stress  will  die  out  in  a power  given  by  the  well- 
known  corner  solutions  referred  to  above. 

A form  of  notch  with  V-shaped  asymptotes,  which  contains  the  parabolic  notch  treated  in  this 
paper,  is  the  following 


r co^(^2)  [“/J/lVs)]  = I_(7r/2)/0  P 


in  polar  coordinates,  where  p is  the  radius  of  the  tip  and  0 = ±3  defines  the  asymptotic  opening  of 
the  notch.  Thus  2 (rr-0)  is  the  asymptotic  angle  of  opening.  The  tip  is  at  a distance 


r = I - I (ir/2 )/3  p 


beyond  the  origin  of  coordinates  on  the  symmetry  axis 
crack  case. 


0 = 0.  If  3 = 7T,  2 r cos-  61 2 = p,  the 


Two  other  points  should  be  mentioned  lest  unfortunate  presumptions  should  be  drawn  from 
the  text  of  this  paper  by  the  reader. 

One,  the  writer  knows  of  no  a priori  reason  to  assume  that  p for  Modes  I and  II  are  equal,  as 
assumed  here  in  the  treatment  for  combined  loading  stresses,  unless  p is  created  to  be  so  artificially, 
as  by  machining. 

And  finally,  the  writer  doubts  that  K calculated  through  use  of  the  formulae  advocated  here, 
will  be  found  to  describe  the  toughness  of  materials  unless  the  materials  are  perfectly  brittle  (-  or 
some  presently  undefined  meaning  or  expression  is  assigned  to  p).  That  is  to  say,  for  example,  that 
if  a fracture  stress  were  calculated  by  the  Mode  I formula  from  the  failure  load  in  a Mode  1 experi- 
ment, the  load  corresponding  to  this  same  fracture  stress,  calculated  by  the  Mode  II  formula  would 
not  be  the  actual  failure  load  in  a Mode  II  experiment,  unless  the  material  was  perfectly  brittle. 


5.  WILLIAMS,  M.  L.  Stress  Singularities  Resulting  from  Various  Boundary  Conditions  in  Angular  Corners  of  Plates  in  Extension. 
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Only  actual  experimental  results  are  likely  to  be  convincing  on  this  point,  but  this  writer’s  approxi- 
mate theory  |2|  and  calculations  indicate  that  the  failures  would  be  sub-surface  and  that  hydrostatic 
stress  associated  with  plastic  deformation,  sufficient  to  bring  the  stress  to  the  failure  level,  is  not 
readily  attained  in  shear.  Nonetheless,  the  complete  theory  of  the  elastic-plastic  interaction  includes 
the  formulae  of  this  paper  as  a limit. 


appkndix 


The  following  equations  are  expressions  for  stresses  at  and  near  ( 1/2  < r/p  < < a/p)  the  tips 


of  crack  like  (a/p  > > 1 ) notches  in  the  directions  of,  and  given  in  parabolic  and  polar  coordinates, 
p is  the  radius  of  the  tip  of  the  nose  and  “a"  is  a length  usually  characteristic  of  its  depth,  being  the 
semi-major  axis  in  the  standard  case  of  an  elliptically  shaped  crack,  while  r is  a polar  coordinate 
such  that  the  tip  is  at  r = p/2,  as  specified  below. 

The  origin  of  the  polar  coordinates  (r,  0)  and  focus  of  the  parabolic  coordinates  as  well  as  the 
origin  of  the  associated  rectangular  coordinates,  is  at  a distance  p/2  behind  the  tip  of  the  nose  of 
the  crack-like  notch,  on  the  crack  axis,  x is  along  the  crack  axis,  the  negative  direction  lying  com- 
pletely within  the  crack.  The  positive  y axis  is  counterclockwise  from  the  positive  x axis.  The 
polar  coordinate  angle  0,  is  positive  in  the  direction  from  positive  x to  positive  y.  The  coordinates 
l(x/a),  (y/a)l , (r/a)  are  oriented  and  centered  in  the  same  manner  as  (x,y),  (r,0)  and  (u,v). 

Polar  Coordinates 

Mode  1 

2 Sr  . 

(p/a)  cos  0/2  . 

- cos3  0/2  + 2 cos  0/2 

(Al) 

kj  Si 

(2  r/a)3/2 

(2  r/a)1/2 

2 S0  _ 
Mi 

h (p/a)  cos  0/2  + 
(2  r/a )3/2 

cos3  0/2 
(2  r/a)1/2 

2Sr0  = 

+ (P/a)  sin  0/2 

cos2  0/2  sin  0/2 

k|  Si 

(2  r/a >3/ 2 

(2  r/a)1/2 

Mode  11 

S[ 

+ (p/a)  sin  0/2 

sin  0/2  (3  cos  0-1 ) 

(A2) 

kllSlI 

(2r/a)3/2 

2 (2  r/a)*/2 

s0 

(p/a)  sin  0/2 

cos  0/2  (3  sin  0) 

k!l  SII 

(2  r/a)372 

2 (2  r/a)1/2 

Sr0  _ . 

f (p/a),  cos  0/2 

cos  0/2  (3  cos  0-1 ) 

kI  Sn 

(2  r/a)3/2 

2(2  r/a)1/2 

Thus  we  have  the  following  for  Sv,  S$  and  Sr0  on  the  boundary  (p/a)  = (2  r/a)  cos-  0/2. 
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Parabolic  Coordinates 


The  crack  boundary  is  a parabola,  designated  by  a parameter  u = uG,  opening  to  the  left.  It  is 
surrounded  by  other  u parabolae.  If  uQ  = 0.  the  sides  of  the  parabola  are  indistinguishable,  though 
separate,  i.e.,  we  have  the  straight  line  crack  (along  the  line  (y/a)  = 0).  Position  points  are  defined 
by  (u.  v)  where  the  values  of  v designate  a set  of  parabola,  orthogonal  to  the  u parabolae,  opening 
to  the  right.  Both  u and  v sets  have  one  and  the  same  common  focus  (origin)  which  is  at  u = 0, 
v = 0.  and  around  which  all  the  parabola  pass,  u2  is  equal  to  the  radius  of  curvature  of  the  tip  of 
the  nose  of  any  u parabola,  as  is  v2  of  any  v parabola.  If  (p/a)  is  such  a radius,  the  tip  is  u2/2  = 
(p/a)/2  to  the  right  of  the  focus  if  the  radius  is  that  of  a u parabola,  and  is  v2/2  = (p/a)/2  to  the 
left  of  the  focus  if  the  radius  is  that  of  a v parabola. 


The  connection  between  (u,  v)  and  [ (x/a),  (y/a)]  is  defined  by  the  mapping  functions  (which 
can  be  incorporated  into  a single  complex  variable  form). 

(y/a)  = u v;  (x/a)  = (u2  - v2)/2  (A5) 

Thus 


u2  = x/a  + r/a  ; v2  = - x/a  + r/a 


( A6) 


where 


r/a  = |[(x/a)2  + (y/a)2]  '/2| 


the  radial  distance  to  the  parabola,  from  the  focus.  From  these  formulae,  geometrical  (as  opposed 
to  algebraical)  construction  of  the  parabolae  is  possible. 

Also. 


i.e., 


(x/a)  1 

" uT  = 1 


['  - 


(x/p)  = J |^1  - (y Ip)2 


u parabolae 


(A7) 


and 


( x/a)  _ 1 
v2  3 


[ Mi  . ,] 


l 


i.e., 


(x/p) 


= | J^(y/p)2  - |] 
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v parabolae 


Evidently,  from  these  expressions,  "a”  as  a unit  of  length  is  not  inconsistent  with  taking  p as 
the  unit  of  length  at  the  nose. 

In  what  follows  we  denote  by  p/a  the  dimensionless  radius  of  the  tip  of  the  boundary  curve 

2 

only,  i.e.,  from  here  on  p/a  = u“. 

The  Mode  1 Airy  stress  function,  Fj,  is  given  by 


2F, 
kj  S| 


u u “ + u-^/3 
o 


(AS) 


and  the  Mode  11  Airy  stress  function,  Fll.  is  given  by 


FI! 

kll  SII 


+ (V  U “ + V u4) 
o 


( A9) 


The  stresses  Su.  perpendicular  to  the  u lines,  Sv,  in  the  direction  of  the  u lines  and  the  shearing 
stress  Suv  are  as  follows. 


Mode  I 


2SU 

(p/a)  u + 

u3 

kl  Si 

u2  + V2)2 

(U2  + V2)2 

:sv 

(p/a) u 

u (u2  + 2 

h si 

(u2  + V2)2 

(u2  + v2): 

2 ®uv 

(p/a) v 

+ V u2 

k|  S| 

(u2  + V2)2 

(u2  + V2)4 

(A  10) 


Mode  II 


(p/a) v 


kII  SH  * (u2  + v2)2  + (u2 


v 


+ v2)2 


Sv 

kn  Sn 


uv 
kII  SII 


— = = T 


+ 0 °/a>  v T 
(u2  + y2)2 

(p/a)  u + 

(U2  + V2)2 


v ( u 2 + 2v2  ) 
(u2  + y2)2 

_ u3 

(u2  + y2)2 


The  top  signs  here  refer  to  positive  0 angles  and  the  lower  ones  to  minus  0 angles. 

Thus  on  the  boundary  of  the  parabola  u = uQ  = |(p/a))/2  I » the  stresses  are  as  follows. 
Mode  I 


2 Su  _ - (p/a)3/2  + (p/a)3 

kI  ^1  (p/a+y2)^  (p/a  + y2)2 


_ + (p/a)3/2  + (p/a)'/2  (p/a  + 2 y2) 

kl  S1  (p/a  + y2)2  (p/a  + v2)2 

_ 2 (p/a)1/2 

(p/a  + v2) 


2 \/^Zp  maximum,  when  v = 0 


2 Syy 

kj  S, 


ifi/jiv + (P/a)  v 

(p/a  + v2)2  (p/a  + v2)2 


(All) 


(A  1 2) 
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' 

''1 


1 


3 ? ? ”1 


J 


. _ 


is. 





Mode  II 


' U = ± (/>/ U ) V T (p/a)  v 

kHSU  (p/a  + v2)2  (p/a  + v2}2 


— 5?L_=  + (phl'L + y (l/a  + 2 y2 

k|*  SH  (p/a  + v2)2  (p/a  + v2)2 


(p/a  + v27 

= + y/^lp.  maximum,  when  v2  = p/a 

— -v—  = + Jil*£!L  + Qp/ap/2 

kH  ^11  (p/a  + v2)2  (p/a  + y2)2 


